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ABSTRACT 


We  examine  a  linear  fractional  flow  model  which  can  be 
interpreted  as  a  Markov  chain  with  partially  controlled 
transition  probabilities.  The  paper  classifies  the  set  L 
of  limiting  distributions  and  details  several  of  its  proper¬ 
ties.  A  precise  classification  in  a  three  dimensional  case 
is  presented. 


1.  INTRODUCTION 


This  paper  identifies  the  set  of  limiting  solutions  for  the  n  dimensional 
constrained  linear  system 

x(t  +  1)  =  x(t)P  +  u (t)  , 

(1) 

x(t)e  =  1  ,  u(t)  _>  0  t  =  0,1,2,  ... 

where  e  is  a  column  vector  with  each  of  its  n  elements  equal  to  one.  The 
initial  vector  x(0)  ^  0  is  given,  and  the  n  *  n  matrix  P  is  nonnegative. 

In  general,  we  assume  w  =  (I  -  P)e  >  0  and  that  (I  -  P)  has  an  inverse. 

Our  principal  result,  however,  requires  slightly  stronger  assumption.  It  fol¬ 
lows  that  x(t)  is  always  a  nonnegative  vector  with  components  summing  to  one; 
i.e.,  x(t)  is  the  distribution  of  some  quantity  at  time  t  .  Equation  (1) 
shows  how  that  distribution  can  change  over  discrete  time. 

Bartholomew,  [1]  and  [2],  has  derived  an  equivalent  expression  of  the 
dynamics  (1)  in  which  x(t)  is  the  distribution  of  a  partially  controllable 
Markov  process.  The  equivalence  is  based  on  the  identity  x(t)w  =  u(t)e  , 
which  holds  if  x(t)  and  u(t)  solve  (1).  For  any  solution  of  (1)  we  define 
z(t)  and  Q[z(t)]  by 

(  x(t)  if  u(t)  =  0 

(2)  z(t)  = 

(  u(t)/u(t)e  otherwise 

and 

Qij(z(t)]  =  Pij  +  w1Zj(t)  , 

or  in  matrix  notation^ 

Q[z(t) ]  =  P  +  wz ( t )  . 

^Since  w  is  a  column  vector,  wz(t)  is  an  n  *  a  matrix. 
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Note  that  z(t)  >  0  ,  z(t)e  =  1  ,  and  that  Q^[z(t)]  is  a  stochastic 
matrix.  It  follows  that 

(3)  x(t  +  1)  =  x(t)Q[z(t) ]  . 

Now,  in  contrast,  suppose  z(t)  is  any  sequence  with  z(t)  ^  0  ,  z(t)e  =  1  . 
Given  x(0)  ,  we  define  Q[z(t)3  and  x(t)  by  (2)  and  (3).  It  is  apparent  that 
u(t)  =  [x(t)w]z(t)  and  x(t)  will  solve  (1). 

This  paper  characterizes  the  set  L  of  limiting  distributions.  We  can  say, 
roughly,  that  for  any  x(t)  and  u(t)  satisfying  (1),  x(t)  converges  geome¬ 
trically  to  L  .  The  set  L  has  two  other  interesting  properties.  First,  let 
a  closed  set  A  be  defined  as  a  trapping  set  if  x(0)  e  A*  x(l)  e  A  .  We 
find  that  L  is  the  smallest  trapping  set.  Second,  if  x(0)  £  L  ,  then  it  is 
not  possible  to  return  to  x(0)  ;  in  contrast,  if  x(0)  is  in  the  relative  in¬ 
terior  of  L  ,  it  is  possible  to  return  to  x(0)  in  a  finite  number  of  periods. 

Section  2  motivates  system  (1)  in  the  context  of  a  manpower  planning  prob¬ 
lem.  Section  3  is  devoted  to  definitions,  a  statement  of  the  theorem,  and  a 
discussion  of  the  result.  In  Section  4  we  examine  a  special  case  with  n  **  3  , 
and  obtain  a  precise  characterization  of  the  set  L  .  Proofs  are  included  in 
Section  5. 

This  paper  extends  and  strengthens  several  results  of  Toole  [7].  Specific 
references  to  Toole's  work  is  included  as  it  appears.  For  completeness  we  have 
included  short  proofs  of  several  of  Toole's  results. 
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2.  MOTIVATION-MANPOWER  FLOW 


Consider  an  organization  with  n  job  classifications  called  ranks.  Let 
be  the  fraction  of  workers  in  rank  i  that  move  to  rank  j  in  one  period 
and  let  v^(t)  1  0  be  the  number  of  workers  hired  into  rank  j  in  period  t  . 
Finally  let  y.(t)  be  the  number  in  rank  j  at  time  t  .  It  follows  that 


ll 

y  (t  +  1)  -  l  y  (t)M  +  v  (t)  , 
J  i=1  *0  J 


or  in  matrix  notation  y(t  +  1)  =  y(t)M  +  v(t)  . 

The  initial  inventory  of  manpower  is  given  by  y(0)  0  .  Assume  the 

organization  is  growing  constantly  at  rate  (0-1)  ;  thus  the  size  at  time 
t  is  0ty(0)e  .  Let  x(t)  be  defined  as  x(t)  =■  y(t)/0ty(O)e  ,  and  define 
u(t)  =  v(t)/0t+1y(O)e  .  Then  x(t)  and  u(t)  obey  Equation  (1)  of  Section  1 
with  P  «  M/6  . 

As  a  second  manpower  example  consider  an  organization  with  n  -  1  ranks. 
Define  y  ,  v  ,  and  as  above.  We  add  rank  n  to  the  organization  to 

consist  of  unfilled  positions,  and  let  y  (t)  denote  the  number  of  unfilled  po¬ 
sitions  at  time  t  .  Define  v^(t)  as  the  number  of  positions  open  during  t 

which  remain  open  in  the  next  period,  set  M  ^  =  0  for  i  =  1,2,  ...,  n  and 
n-1 

-  1  -  £  .  for  i  =  1,2,  ...,  n  -  1  .  With  these  definitions  Equation  (A) 

j=l  ~3 

holds  for  our  second  manpower  system. 
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3.  DEFINITIONS  AND  MAIN  RESULT 


We  define  the  norm  of  a  vector  x  c  ]Rn  to  be  ||x||  =  £  | x . |  The 

j-1  j 

closure,  relative  interior,  and  convex  hull  of  sets  are  denoted  cl  ,  ri  ,  and  C 
respectively.  We  let  <f>  denote  the  empty  set.  The  simplex  S  =  {x  [  xe  =  1  , 
r.  >_  0}  is  the  set  of  possible  distributions  and  we  topologize  the  closed  subsets 
of  the  metric  space  (S,  |  |  ||)  with  the  Hausdorf  metric,  [3]: 


(5) 


6(A,D)  =  Max  Min  | |x  -  y|  | 
xcA  yeD 

d (A,D)  =  Max  [6 (A, D) , 6 (D, A)  ]  . 


Let  E  =  {x  |  x  e  S,x  >_  xP}  be  defined  as  the  set  of  equilibrium  distribu¬ 
tions;  then  the  solution  x(t)  =  x  for  all  t  is  feasible  for  (1)  if  and  only 
if  x  c  E  .  Now  choose  any  z  e  S  and  consider  the  stochastic  matrix 
Q(z)  =  P  +  wz  where  w  =  (I  -  P)e  >  0  .  When  (I  -  P)  has  an  inverse  N  , 
Toole  [7]  has  demonstrated  that  x  -  zN/zNe  is  the  unique  x  r.  S  such  that 

x  =  xQ(z)  .  It  follows  that  x  e  E  if  and  only  if  x  is  a  stationary  vector 

of  some  stochastic  matrix  Q(z)  for  z  e  S  . 

Suppose  x  =  xQ(z)  ,  x  c  S  .  It  does  vot  follow  that  for  any  x(0)  ,  with 

u(t)  =  (x(t)w)z  ,  we  have  x(t)  -*•  x  .  Consider  P  and  z  below. 


(6) 


P 


'0  10 
0  0  1 
h  o  o. 


Z  =  (1,0,0) 


Qtz] 


'0  1 
0  0 
1  0 


=*  X 


(1/3, 1/3, 1/3) 


However  if  x(0)  =  (a,fi,y)  ,  then  x(l)  =  (y,a,6)  ,  x(2)  =  (6,y,a)  ,  x(3)  =  (a,g,y) 
Wien  z  is  strictly  positive,  the  Markov  matrix  Q(z)  is  regular  and  we  then 
have  x(t)  *•  x  for  any  initial  x(0)  . 

Let  A  be  any  nonempty  subset  of  S  <•  nd  R(A)  =  fx  |  x  e  S,x  yP,y  e  A)  . 
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When  A  is  a  singleton,  {>:}  ,  we  use  the  notation  R(x)  .  For  A  t  $  ,  define 

R°(A)  «=  A  ,  R1(A)  =  R(A)  ,  and  for  t  >  l  ,  Rt+1(A)  =  R(Rt(A))  .  For  any  x  e  S  , 

Rt(x)  is  the  set  of  x(t)  feasible  in  (1)  given  that  x(0)  =  x  . 

It  is  easy  to  verify  that  if  A  is  closed,  convex,  or  polyhedral  then  R(A) 

will  have  the  same  property.  Moreover,  A  C  B  =>  R(A)  C  R(B)  ,  and  CR(A)  =  R(CA)  , 

and  cl  R(A)  =  R(cl  A)  .  We  also  can  see  that  for  any  t  ,  E  C  Rt(E)  C  Rt(S)  C  S  . 

Therefore  we  define  the  limiting  set  L  as  L  =  ™  ^(S)  .  Note  that  L  is 

t=0 

nonempty,  closed  and  convex.  Toole  [7]  has  demonstrated 
Proposition  1: 

R(L)  -  L  . 

Proof: 

y  e  R(L)  =*>  y  e  S  and  y  >_  xP  for  some  x  e  L  C  Rt(S)  for  all  t  ^  0  . 

Thus  y  e  Rt+1(S)  for  all  t  0  =*>  y  e  L  .  Conversely  if  x  E  L  ,  then  x  e  S  , 

and  since  x  e  Rt(S)  for  all  t  _>  1  ,  there  exists  a  y(t)  e  Rfc  ^(S)  such  that 

x  _>  y(t)P  .  Let  y  be  an  accumulation  point  of  the  y(t)  .  It  follows  that 

y  e  L  ,  thus  x  >_  yP  =«►  x  e  R(L)  .  ® 

Proposition  2:  (Stanford  [6],  Toole  [7]) 

If  (I  -  P)  has  an  inverse  then 

(i)  int  E  +  *  . 

(ii)  For  any  x(0)  ,  and  y  e  int  E  ,  there  exists  a  finite  T  such  that 

y  e  Rt(x(0))  for  t  >  T  . 

Proof : 

We  have  (I  -  P)  ^  _>  0  .  If  b  >  0  ,  then  y  =  b(I-P)  ^  >  0  ,  and 
x  =  y/ye  e,atisfies  x  e  S  ,  x  >  xP  . 

x ( X  “  p) 

With  x  defined  as  above,  we  have  ?.  -  — — - ,  a  strictly  positive 
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appointment  vector.  Therefore  x(0)Qt(s)  -*•  x  e  ri  E  .  There  exists  an  e 
neighborhood  N  (x)  of  x  such  that  for  any  y  e  N  (x)  we  have  x  e  R(y)  , 

T-1 

or  y  >  xP  .  There  is  a  finite  T  such  that  x(0)Q  (z)  c  N  (x)  .  Thus 

—  e 

T  t 

x  e  R  (x(0))  ,  and  since  x  e  E  ,  we  have  x  e  R  (x(0))  for  all  t  >  T  .■ 

For  k  1  ,  define  the  k*"^  cycle  set  as  =  {x  |  x  e  R^(x)}  .  If 

k  1 

x  e  C  ,  then  it  is  possible  to  return  to  x  in  k  steps;  note  that  C  -  E  . 

% 

co  00  k,  00 

Define  C  =  U  C  ;  if  xeC  ,  there  is  some  finite  k  such  that  it  is 
k=l 

possible  to  return  to  x  in  k  steps. 

Proposition  3;  (Toole  [7]) 

CD  I  00  £  ,  . 

cl  C  =  cl  U  R  (E)  C  L 

Our  theorem  strengthens  this  result  considerably. 

Theorem  1: 

If  w  =  (I  -  P)e  >  0  ,  then 

(i)  L  is  the  unique  closed  nonempty  subset  of  S  satisfying  R(L)  =  L  . 

(ii)  For  any  closed  nonempty  subset  A  of  S  , 

d(Rt(A),L)  £  oCd(A,L) 

where  0  £  a  <  1  . 

(iii)  cl  c“  =  L 

(iv)  If  MACS  ,  then  R(A)  =  A  implies  cl  A  =  L  . 

The  proof  will  be  presented  in  Section  5.  It  depends  in  large  part  on  the 
following  lemma: 

Lemma: 

Let  a  =  [1  -  1/2  Min  (wi  J  i  =  1,2,  n}]  <  1  . 

For  any  x  ,  y  ,  z  e  S  , 
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| j (x  -  y)Q(z) ! |  1  o\ |x  -  y| |  . 

We  suspect  that  a  weaker  version  of  the  theorem  is  valid  under  the  weaker 
hypothesis  that  (I  -  P)  is  nonsingular.  However,  the  lemma  does  not  generalize 
with  that  hypothesis.  Recall  example  (6).  Let  z  =  (1,0,0)  ,  y  =  (0,1,0)  and 
x  =  (0,0,1)  .  We  have 

t  =  0 
t  =  1 
t  -  2 
t  =  3 


This  example  has 

Ik* 

Before  concluding  this  section  we  shall  discuss  several  implications  of 

the  theorem.  First  Rt(A)  ->  I,  for  any  nonempty  subset  A  of  S  .  For  any 

A  we  have  Rt(S)  ^  ^(A)  .  In  addition,  there  exists  a  finite  T  and  x  e  riE 
T  t  t—T 

such  that  x  e  R  (A)  thus  R  (A)  D  R  (x)  for  all  t  _>  T  .  The  sequences 
t  t— T 

R  (S)  and  R  (x)  are  closed  and  converge  geometrically  to  L  .  Moreover, 

t  t-T 

the  sequence  R  (S)  is  contracting,  and  the  sequence  R  (x)  is  expanding. 

As  a  second  point,  let  A  be  any  closed  subset  such  that  R(A)  C  A  .  To 

see  that  L  C  A  ,  assume  x  e  L  and  x  j:  A  .  It  follows  that  x  {  RC(A)  for 

all  t  ,  and  therefore  that 

d[Rt(A),L]  >  d[A,I,l  >  Min  |  |x  -  y|  |  >  0  . 

ycA 


xQt(z) 

yQt(z) 

|  |  (x  -  y)Qt(z) 1 1 

(0,0,1) 

(0,1,0) 

2 

(1,0,0) 

(0,0,1) 

2 

(0,1,0) 

(1,0,0) 

2 

(0,0,1) 

(0,1,0) 

2 

-  y)Qfc(z) ||  =  | | x  —  y | |  for  all  t  . 
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However,  the  term  on  the  left  converges  geometrically  to  zero. 

We  observe  if  x  {  L  ,  then  it  is  impossible  to  return  to  x  .  If  we  did 

k 

return  to  x  in  k  steps  then  x  e  C  C  L  .  To  show  that  we  can  return  to  any 

00 

point  in  ri  L  in  a  finite  number  of  steps  we  must  demonstrate  that  C  is 

oo  Yj 

convex.  If  x  and  y  are  in  C  then  x  e  C  and  y  e  C  for  finite  k 

kli  kh 

and  h  .  This  implies  that  both  x  and  y  are  in  C  .  Since  C  is 

kh  ® 

closed  and  convex,  the  line  segment  joining  x  and  y  is  in  C  C  C  .  It 

follows  Rockafellar  [5],  page  46  that  ri  L  =  ri[cl  C  ]  =  ri  C  C  C  .  For  any 

x  e  ri  L  and  y  e  S  ,  it  is  possible  to  move  from  y  to  x  in  a  finite  number 

of  steps.  This  follows  since  the  sequence  of  closed  sets  RC(y)  -*•  L  ,  and 

T 

x  e  ri  L  ;  thus  there  must  exist  a  finite  T  such  that  x  e  R  (y)  .  In  con¬ 
trast,  if  x  ^  L  and  the  initial  y  e  L  ,  then  it  is  not  possible  to  move  from 
y  to  x  in  a  finite  number  of  steps. 

The  next  section  presents  a  precise  characterization  of  L  in  a  special 
case.  Proofs  of  Theorem  (1)  and  Lemma  (1)  are  presented  in  Section  5. 
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4.  SPECIAL  CASE  - 

■  A  CHARACTERIZATION 

OF  L 

This  section 

examines  the 

special 

.  case 

of 

"PU 

P12 

0  ' 

P  = 

0 

P22 

P23 

0 

0 

P33- 

where  we  assume 

(i) 

wi  > 

0  i  = 

1,2,3 

(7) 

(ii) 

P22  > 

P12 

(iii) 

W2  — 

w^  . 

The  example  corresponds  to  a  three  rank  manpower  hierarchy;  e.g.  assistant,  asso¬ 
ciate,  and  full  professors.  Assumption  (i)  means  it  is  possible  to  leave  from 
any  rank,  (ii)  is  satisfied  if  _>  1/2  for  all  i  and  (iii)  indicates  that 

withdrawal  rates  are  higher  in  rank  2  than  in  rank  3.  We  shall  present  some 
numerical  calculations  below  indicating  that  the  main  result  of  this  section 
is  true  under  more  general  conditions.  The  sets  S  and  E  are  depicted  in 
Figure  1.  In  this  special  case  it  is  possible  to  obtain  a  precise  characteri¬ 
zation  of  the  set  L  . 

For  k  =  1  ,  2  ,  3  let  Q(k)  be  the  matrix  with 


Vk) 


Lk 


+  w. 


P 


ij 


if  j  =  k 

if  J/k 


In  the  manpower  context,  Q ( i )  corresponds  to  making  all  new  appointments  in 

rank  i  .  For  k  =  1  ,  2  ,  3  let  x  be  the  stationary  vectors  of  the  Q(k)  ; 

k  -1 

the  x  are  proportional  to  the  rows  of  (I  -  P)  and  are  the  extreme  points 


of  E  . 
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Recall  that  C  denotes  convex  hull.  Define 

F  =  C{x3Qt(l),x1Qt(2),x2Qt(3)  t  =  0,1,2,  ...}  . 

The  sequence  xiQt(j)  simply  starts  at  x*  and  follows  appointments  in  rank 

2  t 

j  only  for  all  t  .  The  points  x  Q  (3)  for  t  1  all  lie  on  the  line 
2  3 

segment  [x  ,x  )  ,  thus 

F  =  C{x3Qt(l),x1Qt(2),x2  t  =  0,1,2,  ...}  . 


Theorem  2: 

Under  the  assumptions  of  this  section, 

L  =  F  . 

This  theorem  and  the  analysis  developed  in  its  proof  have  two  corollaries. 
Corollary  1; 

For  any  cycle  of  length  k  ,  at  least  one  element  of  the  cycle  lies  in  E  . 
Corollary  2: 

For  any  solution  of  (1)  and  any  c  >  0  neighborhood  of  E  we  have 

x(t)  e  N  (E)  infinitely  often, 
r. 

The  result  allows  us  to  make  an  excellent  and  easily  calculated  approximation 
of  the  set  L  and  to  gauge  the  effect  of  changing  parameters  on  the  set  of 
limiting  possibilities.  Several  cases  are  depicted  below. 

In  Figure  2,  we  have  w^  >  w.^  and  Theorem  2  fails.  However,  in  Figure  3, 
w^  >  W2  and  it  is  obvious  that  our  approximation  is  valid.  Thus  it  is  sufficient 
but  not  necessary  to  have  w^  <  . 

Figures  3  and  4  show  two  alternate  P  matrices  and  the  effect  of  a  change 
in  the  elements  of  P  on  both  the  equilibrium  set  E  and  the  limiting  set  L  . 
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5,  PROOFS 

This  section  details  the  proofs  of  Lemma  1  and  Theorems  1  and  2. 

Proof  of  Lemma  A : 

For  the  moment  consider  z  e  S  fixed.  If  x  =  y  ,  then  |  |x  -  y|  |  =  0  , 
and  the  Lemma  is  trivial.  If  x  /  y  ,  let  ^  for  aH  *  an<*  define 

I+  -  (i  I  >  0}  ,  i"  *-  {i  |  Vi  <  0}  . 

Note  that: 


and 


0 


Ml  -  [v  -[v 

i+  r 


therefore 


IM!  -  2  l  v  -  -2  l  v.  . 

i+  r 


For  each  j  let 


u.  =  )  v,Q..(z)  and  let  J 

J  1 


=  {■ 


uj  >  °)  , 


J  =  {j  |  Uj  <  0}  .  Using  the  same  logic  as  above. 


Mull  -  2  l  u.  =  -2  l  Uj 


In  the  first  case, 
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viQu 


where 

Hence 


If  h+(z)  =  Max  [r, (z)]  , 
ir.I+ 


:  ('  •■)"* 


(z)  «=  h+(z)  |  |  v  |  | 


In  a  similar  fashion,  if 


s.(z)  =  l  Q..(z)  and  h  (z)  "  Max  [s  (z))  , 
J  iel~ 


then 


-2  I  “3  '  -2 


I.vi8i(«)  +  I  V1S1(Z) 

I+  I" 


;  (•• i -,) 


Therefore  | |u|  |  <  h  (z) | |v| |  . 


To  summarize  in  terms  of  x  ,  y  ,  z  , 


!  I  (x  -  y)Q(z)  ||  <  | | x  -  y 1 1  Min 


l  QXj  (z)  .Max  l  Q 

j+  J  r  j" 


For  each 


Qij(z) 


we  have  Q 


ij 


* 

>  w  z . 
3 


where 


Min  {w. 


h  (z)  . 


(z)l 


i  =  1,2,  ...,n}>0 


Thus 
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Min  I  QjjOO  -w  I  zj 
I+  Lj”  J  j" 


and 


Max  I"  I  Qij(z>]  1  1  “  w*  I  zj 

I+  L/  J  j" 


Similarly 


Max  X  Q±j  (z)  <l-w  l  2^  . 

r  Lj"  J  j+ 


It  follows,  that  for  all  z  e  S  , 


Min 


Thus  1 1  (x  -  y)Q(z)  |  [  < 


Proof  of  Theorem  (1); 


Qi1  (z>^ »  Max["  l  Qij(z)| 

i+  Lj+  j  r  Lj-  -* 

(l-^)||x-y||« 


1 1  -  y  <  1 


First  assume 


(8)  d[R(A),K(D)]  <  od(A,D) 

for  all  closed  nonempty  sets  A  and  D  .  It  follows  immediately  that  L  is 
the  unique  fixed  point  of  R  ,  and  with  D  =  L  in  (8) ,  we  obtain  for  each  t  >  1 

d[Rt(A) ,L}  <  otd[A,L]  . 

From  Toole  [7]  we  have  Rfc(E)  C  cl  C  for  all  t  ,  and  C°°  C  L  .  Since  RC(E) 
is  an  expanding  sequence  of  closed  sets  converging  to  L  ,  we  must  have  L  *  cl  C°* 
Finally,  for  any  nonempty  A  with  R(A)  =  A  we  have  cl  R(A)  ■  R(cl  A)  =  cl  A  , 
which  implies  cl  A  =  L  . 

The  proof  is  concluded  by  verifying  Equation  (8). 
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Let  h(u)  =  Min  |  |u  -  v  (  ( 
veR(D) 


Then  6[R(A),R(D)]  -  Max  h(u)  .  Let  u 

ueR(A) 


in  R(A)  be  such  that 


h(u*)  =  6 [R(A) ,R(D) ]  =  Min  ||u*  -  v|j  . 

vcR(D) 


Note  that 


6[R(A),R(D>]  =  6[{u*},R(D) ]  . 

it 

There  exist  some  x  in  A  and  z  c  S  such  that  u  =  xQ(z)  .  Now  choose 

y  e  D  so  that  |  |x  -  y|  |  -  Min)ix  -  y|  |  *»  5[{x},D]  .  Therefore  yQ(z)  e  R(D) 

ucD 


and 


6[R(A)  ,R(D)  ]  »*  Min  |  |u  -v||  <  !  |  (x  -  y)Q(z)  1 1  . 
veR(D) 


From  Lemma  (1), 


6 [R(A) ,R(D) ]  <  o | | x  -  y |  |  =  o5[{x},D]  . 


Also 


5[{x},D]  <  Max  6[{u},D]  *  6[A,D]  . 
ucA 


We  have  shown  that 


5[R(A),R(D)]  <  aS [A,D]  . 


It  follows  in  a  similar  way  that 


6[R(D),R(A>]  <  o6[D,A]  . 


Therefore 


d[R(A),R(D)]  <  od(A,D)  .  M 
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Proof  of  Theorem  (2): 

Since  F  is  nonempty  and  closed,  it  suffices  by  Theorem  1,  to  prove  that 
R(F)  =  F  . 

I 

If  ACS  is  polyhedral  with  extreme  points  u  ,  £  =  1,2,  ...,we  may 
represent  R(A)  as 

C{u*Q(k)  for  £  =  1,2 . k  =  1,2,3}  . 

Therefore 

R(F)  =  C{xV((i  +  l)mod  3)Q(j)  ,  i  =  1,2,3,  j  =  1,2,3,  n  =  0,1,2,  ...}  . 

It  is  clear  that  R(F)  D  F  •  In  demonstrating  R(F)  C  F  ,  our  attendant 
arguments  will  be  sketchy — we  shall  omit  a  large  amount  of  tedious  algebra. 

I.et 

V2  -  {y  :  y  e  S,y2  <  (yP)2> 

and 

V3  =  {y  :  y  e  S,y3  <  (yP>3>  . 

3  t  2  3 

With  the  aid  of  Figure  5,  we  see  that  if  x  Q  (1)  e  V  -  V  ,  we  must  have 

xV(l)Q(3)  e  C{..3,x1Q(3),x3Qt+1(l)} 

and 

xV(l)Q(2)  e  C{x3Q(2),x1Q(2),x3Qt+1(l)}  . 

Similarly,  if  xV(2)  is  in  V3  -  V2  ,  then 

xV(2)Q(3)  e  C{x2Q(3),x1Q(3),x1Qt+1(2)> 


and 
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xV(2)Q(l)  e  C{x2Q(l),x1,x1Qt+1(2)}  . 

Also,  it  is  clear  that 

xV(3)Q(l)  e  C{x3Q(1),x2Q(1)} 

and 

x2Qt(3)Q(2)  e  C{x3,x2}  . 

2 

Now  it  follows  from  the  condition  p22  >_  p^2  that  x  Q(l)  e  F  ,  and  it  is 

easy  to  show  that  we  always  have  x3Q(3)  c  F  •  Hence  we  have  established  that 

3  t 

each  point  of  R(F)  is  in  F  ,  provided  that  the  sequences  {x  Q  (1)}  and 

(x3'Qt(2)}  remain  in  the  sets  V2  -  V3  and  V3  -  V2  respectively. 

It  3  2 

The  condition  p22  £  guarantees  that  if  x  Q  (*.)  e  V  -  V  ,  then 

(2)  xV+1(2)  e  C{xV(2),x2,  (0,1,0)}  C  V3  -  V2  . 

3  t  2  3 

With  the  requirements  p22  _>  p^2  and  w2  >  w^  ,  we  have  that  x  Q  (1)  e  V  -  V 
implies 

(3)  x3Qt+1(l)  e  C{xV(l),x1,  (1,0,0) >  C  V2  -  V3  , 


and  the  proof  is  complete. 


(T‘0‘0) 
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